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ABSTRACT 


Various meteorological parameters were examined to 
determine their predictive value for forecasting intensity 
changes of Western North Pacific tropical cyclones using 
an analog forecast scheme. A number of procedures were 
devised to quantitatively define the forecast potentials 
of these parameters in a systematic method. Based on 
the results of these procedures, predictors for the 24-hour, 
48-hour, and 72-hour forecast periods were defined. 

To evaluate the skill of the predictors and the pro- 
cedures used to define them, an analog intensification 
forecast scheme was developed and applied to a number of 
test cyclones. The results of these forecasts were com- 
pared to the acceptability criteria established by the 
Joint Typhoon Warning Center, Guam and to official fore- 


casts issued by this activity. 
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I. INTRODUCTION 


A. BACKGROUND 

PaoOmeto the mra=.960s, studies in the field of tropi- 
cal meteorology progressed at a very slow pace. A litera- 
ture search through various meteorological journals 
published before 1965 reveals a strikingly small number of 
articles and reports pertaining to tropical meteorology. 
However, since that time, the amount of research in progress 
and the knowledge of the tropical ocean-air regime have in- 
creased by orders of magnitude. Two factors have been pri- 
marily responsible for the increased efforts and results: 
the meteorological satellite and the development of the 
efficient, large capacity, second generation digital conm- 


"se ..ethe advent of the meteoro- 


puter. Garstang [1] states 
logical satellite, and the power of the computer have 
collectively led to great effort and considerable progress 
gm the realm of tropical meteorology...". 

The one ae of tropical meteorology that has been 
the object of more research than any other is the tropical 
cyclone--the hurricane or typhoon. The high winds, large 
storm surges, and excessive amounts of rainfall associated 
with these systems pose a highly destructive threat to many 
densely populated areas of the world, as well as to many 


high density shipping routes and areas of extensive naval 


Operations. Even though tropical cyclones lose their 





organized circulation rapidly as they encounter large land 
masses, the unstable remnants of these systems can cause 
wide-spread damage far inland. Their far-ranging destruc- 
tive nature has been the impetus for extensive research on 
the natural forces which initiate and maintain tropical 
cyclones. 

Much of the research associated with tropical cyclones 
has been directed at improving the accuracy and timeliness 
of tropical cyclone forecasts. The major factors associated 
with the forecasting are the initial development, movement 
Seoeamtensifieatton. With the utilization of high speed 
computers, many objective movement forecast schemes have 
been developed and are itn operational use as forecast aids. 
For example, there are currently six computerized objective 
forecast aids for typhoon movement in use at Fleet Weather 
Central/Joint Typhoon Warning Center Guam’ (JTWC)[2]. None 
of these objective schemes are considered as an adequate 
forecast technique, but are used as aids and are component 
of the final Back, the official forecast. However, there 
are Bees fit 19 no computerized objective schemes in opera- 
tional use as aids to forecasting tropical cyclone intensi- 
fication. The bulk of all forecast research efforts has 
been directed at movement, for if the movement cannot be 
predicted adequately, then even a perfect intensity forecast 
may have limited valve. The overwhelming emphasis placed on 


movement rather than intensity forecasts is documented by 


ee 
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the fact that verification of the latter was not attempted 
Gu IWC prior to 1969. But with the increasing accuracy of 
movement forecasts that has been experienced in recent 
years, the intensity forecast is becoming more relevant [3]. 

Dimes toethe tack of objective aids for forecasting tropi- 
cal cyclone intensification, the forecaster is forced to 
rely on the limited amount of current and recent past data 
at hand to arrive at a forecast. At JTWC the basic forecast 
procedure in use is a linear extrapolation of past intensity 
changes, subjectively modified by expected conditions along 
the forecasted track and by climatology [2]. These proce- 
‘dures are limited in that, for newly developed systems, ex- 
Heapolation is especially unreliable (if at all possible) 
and rapid deepening cannot be forecast to an adequate degree 
of accuracy. An objective scheme ae use as an intensifica- 
tion forecast aid would be a valuable asset to the tropical 
forecaster, 

One of the objective aids for typhoon movement employed 
at JTWC is an Belog scheme called TYFOON. Forecasts are 
provided in the form of probability ellipses of future posi- 
tion,derived from the movement of a number of reports of past 
cyclones that have been screened, selected as analogs, and 
weighted on the basis of their similarity to the subject 
current cyclone. The degree of analog similarity is a func- 
tion of time of occurrence (month and day), geographic posi- 


tion, and a number of meteorological parameters, hereafter 
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called analog elements. The numerical values associated 
with the elements of a particular analog are compared to the 
respective analog elements of the current cyclone to deter- 
mine their resemblance. The degree of resemblance is used 
to assign a weight to the respective analog which will 
define its importance to the final results. A full descrip- 
tion of the original development of TYFOON is given in [4]. 
The overall results of TYFOON have been the best of any ob- 
jective scheme in use at JTWC [2]. The lack of objective 
forecast intensification aids and the success of the TYFOON 
analog scheme for movement were the prime stimulants for 


initiating this study. 


B. OBJECTIVES OF THE STUDY 

The development of the TYFOON forecasting aid necessi- 
tated the selection of various parameters to describe the 
time, position, and synoptic features of current and past 
tropical cyclones in the western North Pacific Ocean. The 
historical data were assembled and compiled into a digital 
format for computer utilization. The primary objective of 
this study was to conduct a pilot study of the analog para- 
meters associated with TYFOON to determine their potential 
as predictors of tropical-cyclone intensification. With the 
assumption that some positive results would become apparent 
from the primary objective, the secondary objective was to 


devise a forecast scheme that would test the predictive 
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value of the analog data and serve as a first attempt at 
providing the forecaster with an objective analog forecast- 


mmcemaid to predict tropical-cyclone intensification. 


meee AREA OF STUDY 

The area of study was restricted to the region for which 
the necessary ae data have been assimilated and digi- 
tized for use by TYFOON. This area consists solely of the 
western North Pacific Ocean. specitically, it is that region 
of the Pacific north of the equator and west of 180 degrees. 
This area experiences more than twice the number of tropical 
cyclone occurrences in comparison with any other area and, 
further, is the only region where the systems appear in all 
12 months of the year [5]. Consequently, the ares is more 
suited to analog forecasting than any other area of tropical 
cyclone occurrence due to the availability of large amounts 


of tropical cyclone data. 


D. DATA 
1. General — 

The data Bah iec for the study were dictated by 
that information which has been compiled for use by TYFOON. 
It consists of various parameters of two major categories: 
cyclone identity elements and analog identity elements. 
These parameters are listed in Appendix A. The data have 
been collected and digitized for each six-hourly report of 
every known tropical cyclone (depression, storm, or typhoon) 


Since 1945, and at the time of this study, such data were 


12 
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available through 1969. The information includes best-track 
data, a result of post-season analysis of operational data. 
The originally observed and analyzed data have been re-ana- 
lyzed to insure consistency in development and movement to 

a degree that may not have been possible at the time of col- 
lection. 

Only limited information was available on the errors 
inherent in the data. Many of the parameters were assigned 
numerical values as a result of suvjective analysis, and 
therefore, it was very difficult to specify the expected 
errors in these data. It is believed that the maximum sus- 
‘tained surface wind has an error of 10 to 15 percent at the 
time of observation and warning, while that error may be 
reduced to 10 percent after post analysis [6]. Estimates of 
errors associated with other data — been eoviaed: and 
are listed in Appendix A. 

When the analog elements for use in TYFOON were 
selected, only cyclone movement was being considered. 
Therefore, the Tee did not include many parameters that are 
Bllgsely associated with typhoon intensity changes, such as 
thermal structure of the cyclone position of the monsoon 
trough, or sea-surface temperature, among others. However, 
there were some parameters such as the past 12-hour change 
in sea-level pressure, the past-movement of the cyclone, 


and the minimum 700-mb height above the vortex, which are 





trom informal discussion with CDR. J. D. Jarrell. 
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known to have some relation to intensity change, at least in 
some situations. 
2. Selection For the Study 

For the purpose of this study, several test cases 
were selected, and for each test case a number of analog 
candidates were defined and compared to the test case with 
respect to intensification characteristics and analog element 
behavior. A test case was defined as one six-hourly report 
of a tropical cyclone. An analog candidate was defined as 
one six-hourly report of any tropical cyclone (other than 
the one from which the test case was selected) which by 
virtue of its position and time of occurrence (month and day), 
was determined to be similar to the test case. 

It is a well dome oc fact that thie rate of occur- 
rence and the general behavior of tropical cyclones vary 
from season ac season and month to month. As to the varia- 
tions in intensification behavior, [7] and [8] provide sub- 
Stantial evidence. Since this study was aimed at the 
relationships between analog Aen nee and intensification, 
it was desirable to eliminate other variables related to 
Mieensification (e.g. month of occurrence) without severely 
restricting the number of possible test cyclones and analog 
candidates. Therefore, it was decided to select for the 
Study only those tropical cyclones occurring in one season 
of the year. It was also considered desirable to initiate 


the research using somewhat well-behaved vorticies. That 
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mma season in which there occurs relatively little rapid 
deepening and little unsuspected weakening over the open 
ocean was considered most favorable for this pilot study. 
However, it was decided that the season selected should 

also contain a relatively large number of occurrences, for 
it has been observed that better results are obtained from 
the TYFOON analog scheme when larger numbers of analogs are 
Meiiable.- Again, it was felt that the most ideal, yet not 
overly restrictive, conditions should be studied initially. 

Studies by Brand [7] indicate that August tropical 
cyclones would best fulfill the above requirements. They 
exhibit relatively little rapid deepening, while also showing 
little weakening over low-latitude open ocean areas. August 
also has the second highest frequency of eeteedee” of 
tropical cyclones. Therefore, it was decided to select at 
least the bulk of the test cases from this month. To allow 
for a variable time span from which analog candidates co ld 
be selected for any August test case, it was necessary t 
mmelude the tropical cyclones of July and September in tie 
data. 

Bacon on the above considerations, the analog data 
for the months of July, August and September for the years 
1945 to 1969 were obtained (on magnetic tape) from the 
Environmental Prediction Research Facility, Monterey, Calif- 


Ornia. 


Proomeinrormal diseussion with CDR. J. D. Jarrell. 


Stealer of six-hourly reports. 
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ii. PROCEDURES 


Peeper INTLTEON OF CYCLONE INTENSLTY 

In a study of tropical cyclone intensification, the 
first necessary step is to define intensity. The variables 
most commonly associated with the various stages of develop- 
ment of a tropical cyclone from depression to storm to 
typhoon or hurricane include the amount, orientation, and 
maximum height of clouds, the minimum sea-level pressure, 
and the maximum wind. The content of the available data 
‘eliminated the cloud conditions as the intensity parameter 
for this study. Of the remaining two parameters, minimum 
sea-level pressure was estimated to be subject to smaller 
observational errors than maximum wind, and thus would be 
a more reliable indicator of intensity. However, it is the 
maximum wind that is of prime interest to the forecaster 
and recipients of his services; therefore, it was this para- 
meter that was selected to describe the tropical cyclone 


intensity. 


Be. SELECTION OF TEST CYCLONES 

For the initial study, a total of 30 test cases were 
selected. The selection was not random, but was carried 
out in such a manner to insure that a wide range of initial 
intensities and positions were represented. Widely varying 


intensification characteristics, to include development and 
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weakening, were also exhibited by the test cases. As men- 
tioned before, it was not desirable to be overly restrictive 
with the data as this would possibly inhibit favorable re- 
sults. Based on the findings in this study, restrictions 
could be applied, if necessary, in subsequent studies. In 
other words, if one adequate objective forecast scheme could 
not be derived for cyclones of widely ranging initial in- 
tensities and positions, then the input data could be strati- 
fied with respect to these variables and different forecast 
schemes attempted for the various classes of cyclones. 

Other factors were considered in the test cyclone selec- 
tion. Since one avenue of research utilized historical data 
of the test cases, a restriction was placed on the minimum 
amount of history possessed by a number of the cases. There 
were, however, some test cases with no recorded historical 
data. Recurvature was also considered and a number of the 
test cyclones were selected on the basis of whether or not 
recurvature took place within the forecast periods or in the 
Becent past. 

As stated ee” ene bulk of the test cases were 
selected from the month of August. However, to insure that 
a sufficient number of appropriate test cases from recent 
years collectively represented a wide range of initial in- 
tensities and positions, and to avoid an excessive number 
of related test cases (from the same tropical cyclone) it 


was found necessary to select a few test cyclones from the 
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month of September. The 30 test cases initially selected 
Memeethe study are listed in Appendix B. 

The total number of six-hourly reports included in the 
three-month period of data obtained for the study amounted 
to over 7500. It is realized that only 30 oe cases may 
not constitute a representative sample of all analog para- 
meters involved. However, the full range of possible 
initial intensities and a highly variable distribution of 
initial positions were represented. It was believed that 
if useful persistent relationships (between the analog 
elements and intensification) not only exist, but also re- 
main valid for highly oe initial conditions, then at 
least some indications should be observed from as few as 
30 test cases. Again, it is emphasized that the prime 
purpose of this study was " Cetermane ee some, analog 
elements possess any potential as predictors of intensifi- 
cation, not to quantitatively describe the various degrees 


of this potential. 


C. BASIC APPROACH TO THE MAIN OBJECTIVE 

The main objective of this study required the develop- 
ment of a procedure that would examine the analog elements 
and identify any behavior which appeared to be related to 
intensification. Several trial and error retest schemes 
could have been attempted using various parameters with 
associated weight factors being assigned intuitively. It 
was decided, however, that a systematic approach would be 


more likely to bear fruitful results. 
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In an analog intensification forecast scheme, a compari- 
son ‘would be made of known elements (predictors) of the 
current cyclone with those of past cyclones that occurred 
within specified space and time limits, as well as exhibiting 
an intensity within some predetermined range. Based on a 
comparison of the analog elements, a number of analog candi- 
dates that most closely resemble the current cyclone would 
be determined. The changes in intensity exhibited by these 
most similar analogs would be formulated into a forecast for 
the current cyclone. 

To investigate the forecast potential of the analog ele- 
ments, it was decided to use the inverse of an analog fore- 
cast procedure. For a given test case of known future 
intensification, all six~hourly reports of past cyclones 
were screened to determine those that occurred within speci- 
mared time, space, and initial intensity limits. Then, 
rather than classify the analogs on the basis of their 
element values (as would be done in a forecast scheme) the 
future intensification of each selected analog candidate was 
Meinined, and anes were classified according to their 
deviations in future intensification from that of the cur- 
rent cyclone. A class was defined as a collection of analogs 
that exhibited similar deviations in intensification from 
that of the test case for a given forecast period. Each 
class of analogs represented an intensification range, and 
the similarity of that range to the test cyclone intensifi- 


cation varied from class to class. The analog elements 
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associated with the candidates in each class were ranked to 
determine their average degree of matching with respect to 
the test case elements. This basic approach will become 
better understood as the various procedures are discussed. 

The hypothesis around which the basic approach te the 
study was designed specifies that if any of the analog ele- 
ments were related to intensification, then the degree of 
Matching of these elements, with respect to the current 
cyclone values, should decrease as the deviation in inten- 
sification of the associated analog candidates, from that 
of the test case, increases. If the deviations in inten- 
sification are thought of as forecast errors, then the 
hypothesis implies that higher degrees of matching are 
associated with smaller forecast errors for potential pre- 
dictors of intensification. If any of the elements 
exhibited this trend with some degree of consistency, then 
a class of analogs with intensification similar to the test 
case could be determined, and this intensification used as 
a forecast aid. This basic procedure was considered to 
be the most erie vent and direct approach to the main objec- 
tive. 

The main problems encountered in this basic approach 
were deciding on what criteria to use for the analog classi- 
fication and what technique to use for the degree of matching 
computation. The following sections describe the various 
specific procedures that were developed in an attempt to 


accomplish the primary and secondary objectives. All 
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procedures were programmed in FORTRAN IV for use with the 
IBM 360 computer at the Naval Postgraduate School, Monterey, 


@arifornia. 


Dee EF ERST APPROACH TO THE OBJECTIVE 

The first procedure developed for the study was an in- 
verse analog program, as described in the previous section. 
It consisted of two basic operations: (1) test case defini- 
tion and analog selection; (2) analog classification and 
parameter rank computation. 

In the first operation a test cyclone was defined by 
means of a test case identification number, by its associ- 
ated latitude, longitude; and month, day and time of occur- 
rence, and by its associated analog elements. The future 
miimemsifications for 12, 24, 48, and 72 hours, exhibited 
by the test case, were also defined. For test cases not 
Seep iting a particular future intensification, a code of 
999 was assigned. The intensification for each forecast 
period was defined as the difference between the intensit 
at the end of the respective period and the initial 
me nsity. 

The next step was to define the time, space, and initial 
intensity limits to be utilized in the analog selection. 
Time and space screening criteria of plus or minus 15 days 
from the day of the test cyclone, and pius or minus five 
degrees latitude and longitude from the initial position of 


the test case were selected for this procedure. In general, 
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vorticies of extremely different intensities exhibit unlike 
intensification characteristics over an extended period of 
time. For example, a young developing tropical depression 
would be unlikely to exhibit intensification characteristics 
similar to that of a fully-developed super typhoon. Under 
conditions favorable for development, the young cyclone 
could continue to intensify for a considerable length of 
time, whereas the super typhoon, being near or at its point 
of maximum intensity, would not be expected to exhibit this 
behavior. In an analog forecast scheme, an initial intensity 
screen would be utilized to eliminate unnecessary and un- 
‘desirable noise in the selected analogs. Since the analog 
selection phase of this procedure was similar to that which 
would be utilized in an analog forecast scheme, the initial 
intensity screen was applied. Limits of plus or minus 20 
knots from that of the test case were defined for the analog 
selection. With the test ee and its associated screen- 
ing criteria established, the procedure continued with the 
analog Pelacesern, | 

Each sieeneeerly ae of all past tropical cyclones 
(contained in the data), excluding the cyclone with which 
the test case was associated, was screened with respect to 
its time of occurrence (month and day), position, and initial 
intensity. All those reports satisfying the previously 
defined screening criteria were designated analog candidates. 
For each analog selected, the associated 12-, 24-, 48-, and 


72-hour future intensifications were determined. With the 
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analogs selected and their intensifications computed, the 
first operation in this procedure was completed. 

The analog classification and parameter rank computation 
proved to be the most difficult aspect of the study. No 
previous guidelines were available to indicate what would be 
the best approach. The various methods of attack were some- 
what trial and error, but all were based on the hypothesis 
Stated earlier, and all had one concept in common, namely 
the analogs were classified according to their deviations 
in intensification from that of the test case for each fore- 
cast interval. 

In the first approach to the analog classification, 
positive and negative deviations ware considered. For each 
forecast period the following analog classification scheme 


was implemented: 


GCiass 1 Deviation = O>- +9 KTS 
Class 2 Deviation = O>- -9 KTS 
Class 3 Deviation = 10> 19 KTS 
Class 4 Deviation = -10->-19 KTS 
Class 13 Deviation = 70> 79 KTS 
Class 14 Deviation = -/0--79 KTS 


With all (or most) of the analog candidates classified, the 
next step was to compute analog element ranks. Analogs 
with deviations in intensification not defined above were 
eliminated from further investigation. 

For each analog element associated with an analog can- 
didate in the above classes, a rank was computed to deter- 


mine the degree of matching with the respective test cyclone 
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element. In this approach the ranges of the analog elements 
Over the entire data file were Be eee ee The rank of the 
hl 


N element of each classified analog candidate was computed 


using the following formula: 


PA(N) - PT(N) 


RCN) = (1.0 ap geile Ova; (1) 
eee 

PA(N) = Noh element value of the analog 

mi{(N) = N element value of test cyclone 

r(N) = range of Nth element 

R(N) = rank of Nth element associated with the 


particular analog 


PA(N) - PT(N) 


The ratio r(N) 


indicates how different the two 
values of the iad element are as compared to the maximum 
Meosaible difference. Subtracting this quantity from 1.0 
results in a higher rank for a closer match. For an identi- 
Seaeematch, PA(N) = PT(N) and R(N) = 100.0: a rank of zero 
would indicate the two values, PA(N) and PT(N), were on 
Opposite extremes of the range. 

Once the elements were ranked for all classified analo s, 
the average rank of each element was computed for each clé:3s, 
Bus arriving at an average degree of matching. The standard 
deviation of the individual ranks in each class was also 
calculated for each analog element. Also for each element, 


a graph of average rank versus class deviation in intensifi- 


cation was plotted. Recalling that these deviations could 


>Max imum and minimum values of the analoe elements and the 
ranges were computed in a separate program and specified 
in this procedure. For direction of movement elements, 
the range was specified as 180 degrees. 
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be considered forecast errors, these graphs were essentially 
the average degree of matching versus forecast error. 

This procedure was applied to the first 15 test cyclones 
listed in Appendix B, and the results were analyzed. Graphs 
revealing a decreasing average degree of matching with in- 
creasing forecast error were not observed with any persis-— 
tence for any of the analog elements. In several cases, an 
increasing average degree of matching with increasing fore- 
cast error resulted with the past direction and speed of 
Movement elements. No physical interpretation could be 
given for this behavior. 

One obvious fault in the procedure was noted. For the 
48- and 72-hour forecast periods, and to a lesser extent the 
24-hour period, the selected analogs in any one class ex- 
hibited markedly different intermediate intensification 
characteristics. In a 48-hour period, for example, one ana- 
log may have intensified 30 knots in a gradual but steady 
manner, while another analog may have intensified 50 knots 
in the first 24 hours, remained steady for 12 hours and then 
Peeicened by 20 knots in the next 12 hours. These two hypo- 
thetical cases exhibit markedly different behavior, but 
Since each experiences a net 48-hour intensification of plus 
30 knots they would be classified together under the current 
procedures. Since it was the intention to classify only 
Similar analogs together in order to study the variable be- 


havior of the analog elements, modified definitions of 
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intensification were implemented for the 24-, 48-, and 72- 
hour forecast periods, and are shown below. The letter I 
represents intensity and the subscripts indicate the number 
of hours from the time of occurrence of the test cyclone 


or analog. 


4 15, = 194 7- Fre 
S Eggs pam yy (2) 
Bie 


For purposes of analog classification for 24-hour intensi- 
fication, the deviation in intensification from that of the 
test case between time = 12 and time = 24 had to fall within 
Peeoet of class limits. For 48-hour intensification, the 
deviations for each 12-hour period involved (24 to 36 and 36 
memaeeso hours) Had to fall within one particular set of class 
limits. Similar restrictions were placed on the analog 
classification for the 72-hour intensification. Classifica- 
tion for the 12-hour period remained unchanged. Compared 
meomtne Original definitions of intensification, these defini- 
tions would require more time-restricted relationships to 
exist between the analog elements and intensification in 
order for the desired trend of decreasing degree of element 
Matching with increasing forecast error to appear with any 
persistence. As an example, for an element to exhibit a 
Semone relationship to 48-hour intensification, the relation- 


ship would have to be valid during more precisely defined 
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mime Limits (24 to 48 hours). Under the prior definitions 
@emancensification, an indicated relationship between the 
behavior of an element and the 48-hour intensification 

could have been actually very strong for the first 24 hours, 
but relatively weak for the second 24 hours. However,these 
restrictions were considered more favorable than to contend 
with the identical classification of analogs exhibiting 
highly variable intensifications. 

The first 15 test cases were rerun using the modified 
definitions of intensification, and the results evaluated 
again. Little if any change from the previous results were 
observed. In most Be oi the average ranks of all the 
elements were very high, generally greater than 85. This 
was due to the use of the element ranges in the equation 
for the ranks. The computed ranges were usually the result 
of extreme and uncommon element values. For example, the 
range of sea-level pressure extended from 875 mb to 1012 r .. 
To compensate for this, the standard deviations of the 
analog element distributions over the entire data file we e 
computed and substituted for the ranges. With the use of 
the standard deviations, negative ranks would result when 
the difference in values of an element was greater than the 
standard deviation of the element distribution. A closer 
match of an analog element to its test cyclone value would 
result in a more positive rank. 

A third run of the first 15 test cases revealed no 


positive results. Plots of the average rank versus forecast 
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error revealed a much more variable pattern than previously 
observed, but no consistent decreasing degree of matching 
with increasing forecast error resulted for any of the analog 
elements. Many cases revealed extreme oscillations of the 
average ranks of some elements from one class of analogs to 
the next. 

Only one relationship persisted throughout the entire 
Study to this point. In almost every instance for any analog 
element in any analog class, a larger average rank was 
associated with a smaller standard deviation of the element 
ranks involved. Only a small increase in an average rank 
was generally accompanied by a small decrease in the standard 
deviation. Therefore, lower average ranks were associated 
with relatively more variable individual ranks, but higher 
average ranks were associated ei fh variable individual 
ranks. 

The fact remained, however, that both high and low 
average ranks were observed in all classes of analogs for 
all 15 eC It was felt that if some persistent and 
useful trend were to develop, at least some weak indications 
would have resulted by this point in the research. In the 
absence of any such indications, this method of analog 
classification and element ranking was abandoned. The possi- 
bility existed that the ranking procedure utilized was too 
sensitive for the accuracy of the data involved and for the 


reliability and persistence of some relatively weak, but 
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somewhat useful relationships which might exist. It was 
decided to attempt a less restrictive classification pro- 


cedure and to adopt a less refined ranking technique. 


E. SECOND APPROACH TO THE OBJECTIVE 

The two basic operations of the first approach were 
also utilized in this procedure. The test cyclone defini- 
tion and analog selection phase was modified in two ways. 
Pmerst, the initial intensity limits for the analog screen- 
ing were expanded to plus or minus 40 knots from the 
initial intensity of the test case. It was felt that a 
greater number of analogs would generate better results, as 
has been observed with the TYFOON analog scheme. Due to 
possible diurnal effects on the results, the 12-hour fore- 
cast was excluded and the 24-hour intensification re- 
defined as 


A To, 7 ta, 7 Lo (3) 


The analog classification and element rank computation 
were completely re-developed for this procedure. For analog 
@iasesification, the deviations in intensification during 
the 24-hour periods or ictal times pomeoehours, 24 to 43 
hours, and 48 to 72 hours were considered. The deviations 
within the individual 12-hour periods contained in these 24- 
hour periods were not considered in order to make the classi- 
fication less restrictive and to eliminate possible diurnal 


effects. For each forecast period (24-, 48-, and 72-hourS), 
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the analogs that intensified within plus or minus 10 knots 
of the test case per 24-hour period were identified and 
classified together. 

The next step was to assign a rank to each element of 
the classified analogs. Rather than define a degree of 
matching for each individual a emenen a technique was de- 
veloped that would indicate a match or non-match result 
for each element and all matches would be given equal 
weight. It was necessary to establish for each analog 
element an interval, centered on the test case value, within 
which an analog candidate value would have to occur to be 
considered a match. Lt a match occurred, the element of 
Bile analog was assigned a rank of 1.0. If a match failed 
to occur, the element was assigned a rank of 0.0. 

Due to the variable ranges of possible values among 
the analog elements each established interval of acceptance 
had to be unique to a particular element. Since the 
Standard deviation of a data population indicates the va ia- 
Peetity of the quantity, the interval of acceptance was 1m ide 
Proportional to siti measure of dispersion. An interval 
of plus or minus one fourth of a standard deviation centered 
on the test case value was selected as the interval of 
acceptance for each analog element. All analog element 
values occurring within this range were considered a match; 
all others were rejected. The acceptance intervals for the 


analog elements are given in Appendix C. 
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For each element the number of occurrences of missing 
data associated with the classified analogs was determined. 
Subtracting these numbers from the total number of analogs 
in the class resulted in the number of possible matches 
for each element. By summing up the ranks associated with 
a particular element over all the classified analogs, the 
number of matches was determined. Dividing this by the 
number of possible matches and multiplying by 100 resulted 
in the percent number of matches for the element. This 
quantity was now referred to as the degree of matching. 
This procedure was applied and a degree of matching estab- 
lished for each of the 18 elements. 

The classification procedure was now repeated for all 
the analogs that exhibited an intensification within plus 
or minus 20 knots of the test cyclone per 24-hour period, 
and for each forecast period. The degree of matching of 
each element was computed exactly as before. The procedures 
were again repeated for deviations within plus or minus 30 
knots and finally for deviations within plus or minus 40 
knots per DTN. period, for each forecast period. 

For each element, the degree of matching versus the 
Maximum allowable class deviation in intensification (maxi- 
mum allowable forecast error) was plotted. As in the first 
approach a decreasing degree of matching with increasing 
forecast error was considered to be indicative of a relation- 
ship between element behavior and intensification. 

As the maximum allowable forecast error increased from 


10 to 40 knots the number of analogs and therefore the number 
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of possible matches inereased. But for a potential predic- 
tor, any increased number of actual matches was expected to 
reSult in a lower degree of matching. 

The procedure was applied to all 30 test cases listed 
in Appendix B, and the results analyzed. The desired trend 
of decreasing degree of matching with increasing forecast 
error was exhibited to at least some extent by all analog 
elements but much more frequently by some elements. 

The next step was to determine the relative signifi- 
cance of each analog element as a potential predictor of 
intensification. For each forecast period and for each of 
the analog elements, the number of occurrences of decreasing 
degree Ofmatch with increasing forecast error was noted and 
the percent of this number to the total possible number of 
occurrences was computed. These percentages were referred 
to as indicated forecast potentials. The total possible 
number of occurrences varied with individual test cases. 
some cases contained missing data which ruled out possible 
matches for some elements and in some cases a future 48- 
or /2-hour verifying he enetee Gia mot e€xist whieh prohibited 
analog classification. The results of these computations 
are presented in Appendix D. 

As a first attempt, it was decided to choose for each 
forecast period the analog elements associated with the 
five largest indicated forecast potentials as predictors 


of intensification for the respective period. 
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Before developing an experimental forecast scheme, one 
further step was conducted to provide a more rigorous study 


of the analog elements. 


feet HiRD APPROACH TO THE OBJECTIVE 

In order to make the study as thorough as time would 
allow, various parameters were combined to derive other ana- 
log elements to be examined for forecast potential. Also, 
since it was the change in intensity that was to be forecast, 
the recent time rates of change of some original and derived 
elements were investigated for forecast potential. Only the 
past 24- and 48-hour changes of various parameters were 
examined, as utilizing any greater amount of history might 
lead to limitations in a forecast scheme applied to a re- 
cently detected cyclone. The derived elements are listed in 
Appendix E. All but three were easily derivable from the 
original analog and identify elements of TYFOON. 

The three new elements which required some external in 
formation were the’ 1000/700-mb thickness elements. These 
elements were included as indicators of the thermal charac- 
teristics of the lower domain of the cyclones. The minimum 
700-mb height above the cyclones and the minimum sea-level 
pressure (original analog elements) made possible an approxi- 
mation of the 1000/700-mb thickness. 

The missing term for the thickness elements was the 
height of the 1000-mb surface. It is a widely-held belief 
that a sea-surface temperature of at least 26-27C is required 


for tropical cyclone development [5]. The air temperature 
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memar the air-sea interface musSt be similar to the sea-sur- 


Face temperature. 


A mean virtual temperature of 30C for 


the layer between the surface and 1000 mb was considered a 


reasonable average approximation. With this approximation, 


and the use of the hydrostatic equation and the equation of 


State, the height of the 1000-mb surface could be estimated. 


OP = -gpdz (4) 
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The 1000-mb height thus obtained was subtracted from the 


minimum 700-mb height above the vortex to determine the 
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thickness. In this manner, the current, past 24-hour and 
past 48-hour thicknesses were computed and the appropriate 
time rates of change calculated. 

The investigative procedures conducted on these derived 
elements were essentially identical to those applied to the 
Original elements and described in the previous section. 
Some program modifications were necessary, however. 

In the test cyclone definition and analog selection 
phase, the necessary past 24~ and 48-hour data were defined 
and a routine was added to compute the time rates of change. 
A history indicator was also defined for each test cyclone 
to indicate the extent to which history elements of the 
selected analogs should be calculated. This phase was also 
modified to compute the necessary and available derived 
elements for each selected analog. 

Only one modification was necessary for the analog 
classification and element rank computation phase. It was 
necessary to redefine the matching criteria (acceptance 
intervals) for the new elements. To compute the standard 
Mee ations of all the elements would have required more com- 
puter time than was practically available. Acceptance inter- 
vals were defined subjectively based on the estimated 
variability of each element and on the estimated probable 
errors in the values. As before the intervals were centered 
on the test cyclone element values, and for an analog 


~ 


element to be accepted as a match, its value had to occur in 
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mie appropriate interval. The acceptance intervals for the 
derived elements are listed in Appendix E. 

The procedures were carried out on the 30 test cases 
listed in Appendix B and an analysis of the results was made 
as described in the previous section. For each forecast 
period, the indicated forecast potential was computed for 
each element, and the elements associated with the five lar- 
gest indicated forecast potentials were selected as predic- 
tors of intensity. The indicated forecast potentials of 
the derived analog elements are listed in Appendix F. 

WeEietne tntenscwty predietors defined, it was now neces-— 
Sary to establish their relative importance for use ina 


forecast scheme. 


fee DERIVATION OF WEIGHT FACTORS FOR SELECTED PREDICTORS 

The results of the study to this point revealed that the 
potential among the original and derived analog elements as 
Meecaictors of tropical cyclone intensification varied 
greatly. For the elements that were selected as predict srs, 
in a forecast scheme, more emphasis would be placed on those 
that exhibited more persistently the trend of decreasing de- 
gree of matching with increasing forecast error. The degree 
of emphasis placed on each predictor would be indicated by 
a weight factor. Sit sol wGen predayctors, for example, it 
would not suffice to assign the best one a weight of ten and 
the worst one a weight of one. This would imply that the 


best predictor possesses ten times the forecast potential of 


36 








the worst one, whereas it may only possess twice as much. 
Therefore it was necessary to derive weight factors as a 
function of the indicated forecast potentials resulting 
from the earlier investigations. 

To assign weight factors relative to the indicated 
results, the following procedures were applied. First of 
all, a small number of selected predictors were eliminated 
due to their close relationships to other predictors with 
equal or higher indicated forecast potentials. This avoided 
a concentrated emphasis on the related predictors, Eee Wl 
lative effect of which, might often lead to poor results. 
For each forecast period the selected predictor (original 
or derived) with the lowest indicated forecast potential 
was assigned a weight factor of 1.0. The weight factors 
for the other predictors for each period were computed by 
dividing the values of the associated indicated forecast 
potentials by the lowest potential for the respective fore- 
cast period. This procedure normalized the numerical weight 
factors based on the smallest value of an indicated forecast 
potential for each period. The selected predictors and 
their meeipective weight factors are listed in Appendix G. 

At this point, the primary objective of the study was 
accomplished. An examination of the potential value of the 
analog elements as predictors of intensification was made; 

a number of predictors was selected for each forecast period 
and associated weight factors derived. To test the findings 


of this part of study and to evaluate the merit of the 
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procedures utilized for identifying useful predictors, a 


forecast scheme was developed. 


H. DEVELOPMENT OF A FORECAST SCHEME 

The forecast scheme to be developed could function in 
many different ways. A weighted average intensification 
of all the selected analogs could be computed as a forecast, 
or a filtering technique could be applied to all the analogs 
to determine those which most closely resembled the test 
cyclone and apply the intensifications of these filtered 
analogs as a forecast. A scheme of the latter type was de- 
vised, as the actual forecast potentials of the predictors 
memed be more directly utilized. The following forecast 
scheme was developed and applied to the 30 dependent cases 
used to select the predictors. 

The test. cyclone definition and analog selection opera- 
tion was identical to that developed for analyzing the de- 
rived elements with one addition. For each forecast period 
(0-24, 24-48, 48-72 hours), the maximum positive and nega- 
tive intensifications exhibited by the selected analog can- 
didates were determined. 

The second operation constituted the actual forecast 
scheme. The 24-hour forecast was first determined. All 
selected analogs were classified according to their 24-hour 
intensification. The analogs that intensified between —-5 
and +5 knots, inclusive, were classed together and assigned 


emeiaracteristic intensification of zero. Next, all analogs 
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that intensified between +5 and +15 knots were classed to- 
gether and assigned a characteristic intensification of 
foeeknots. Analogs that exhibited an intensification of 

+5 knots were assigned to both classes as the data were 

not accurate enough to warrant a rigid delineation between 
classes. The procedure was repeated until all the inten- 
sifying analogs were classified. Next, all the analogs 
exhibiting negative intensification were classified in a 
Similar manner. One of the resulting characteristic inten- 
sifications would be selected as a forecast. 

The characteristic intensifications and the number of 
analogs classed with each were printed out. The distribu- 
tion of the analogs cee the classes indicated the most 
bo Least frequently occurring intensifications of all past 
Ssix-hourly reports of all past tropical cyclones exhibiting 
a time of occurrence, position, and initial intensity sim?- 
lar to those of the test cyclone. The epee eon there - 
fore represented an intensification climatology or a 
precisely ae eueiea group of tropical cyclone reports. The 
Characteristic intensification of the class containing the 
erect number of analogs was defined as the modal intensi- 
mreation. 

For each predictor, the degree of matching was deter- 
mined for each analog class. As before, the degree of match- 
ing was defined as the number of actual matches relative to 


the number of possible matches. The matching criteria 
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(acceptance intervals) were those utilized in the selection 
of the predictors. The analog class that exhibited the 
highest degree of matching for each predictor was determined, 
and the weight factor associated with that predictor was 
assigned to that class of analogs. If two or more classes 
exhibited the same degree of matching for a predictor, the 
weight factor was assigned to all of those classes. The 
assigned weight factors for each class were then summed and 
the summation for each class was printed out. These same 
precedures were then applied to the 48~ and 72-hour forecast 
periods. 

Originally, it was the intention to select the charac- 
teristic intensification associated with the largest summa- 
tion of weight factors as the forecast. However, analysis 
of early results revealed a apenian area and a previously 
Meese ooked facter: the variable probability of chance ele- 
ment matching. The numbers of analogs associated with the 
Meeererent classes varied from as few as one to ten,to as 
Many as 80 to 100, as an example. In most cases, the 
misses containing relatively very few analogs resulted in 
relatively very large weight factor summations. The trend 
was so persistent that some unaccounted-for factor was sus- 
pected to be influencing the results. It was believed 
that the smaller the number of analogs in a given class, 
the greater was the probability of getting element matching 


due to pure chance. The implications of this factor are 
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discussed further in the results section and only the pro- 
cedures attempted to account for chance matching are stated 
here. 

To compensate for the chance matching factor, it was 
necessary to account for the number of analogs in each class. 
Simply multiplying the weight factor summation of each class 
by the number of associated analogs would not suffice, be- 
Cause the ranges and orders of magnitude of the numbers of 
analogs were much larger than those of the weight factor 
summations. The results would be almost totally functions 
of modal intensification. To reduce the effects of the num- 
ber of analogs in each class on the results, the square 
roots of these numbers were computed and multiplied by the 
respective weight factor summations. This procedure alle- 
viated the problem of chance matching to a large degree, 
but in many cases in which there were extremely large varia- 
tions in the numbers of analogs among the classes, the prob- 
lem still persisted. Therefore it was decided to eliminate, 
as possible forecasts, the characteristic intensifications 
of classes with relatively very few analogs. In order to 
Maintain a strictly objective forecast scheme, a definite 
procedure was needed to select those classes whose charac- 
teristic intensifications would be retained for the forecast 
determination. Only those characteristic intensifications 
whose frequencies of occurrence were in the upper 50 percent 
Of the entire distribution were retained. That is, the 


classes were selected in order of decreasing number of 
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associated analogs until one half of the total number was 
selected, unless an equal number of analogs in two or more 
classes required more than one half of the total. The 
development of this procedure was purely subjective, but 
it did provide for an objective determination of a forecast, 
and a complete revision of the forecast scheme in order to 
account for chance matching was impractical at this time. 
For each test case, three forecasts were derived based 
on different criteria: 
Duetiiodwl meen siricati1on: characteristic 
intensification associated with the 
largest number of analogs 
Pieewelentred predictors: characteristic in- 
tensification associated with the largest 
weight-factor summation (after elimina- 
tion of specified analog classes) 
Epeewetlented climatology: characteristic 
intensification associated with the 
largest product of weight-factor summa- 
tion and square root of the number of 
associated analogs (after elimination 
of specified analog classes). 
For each method the 24-hour intensification was determinec 
and added to the initial intensity to obtain the 24-hour 
forecast intensity. Then the 24- to 48-hour intensification 
was determined and added to the 24-hour intensity to obtain 
the 48-hour forecast intensity. The 72-hour intensity was 
determined similarly. A minimum forecast intensity of 30 
knots was imposed on each forecast method. This was done to 


reduce the degrading effects of a highly underforecasted 24- 


hour intensity, for example, on the 48- and 72-hour forecast 
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intensities which might otherwise be very accurate. The 


resulting forecasts were compared to the official JTWC 


forecasts issued for the 30 
was then applied to 62 more 
uSing the same criteria for 


These test cases are listed 


test cases. The forecast scheme 
independent test cases selected 
the original 30 test cases. 


in Appendix H. The results of 


these forecasts were also compared to official JTWC fore- 


casts. DMOGeuSSiton Of EMe £Orecast results and the results 


Othe Study in general is given in the next section, 








aie ReEoULTS AND CONCLUSIONS 


PDL SCUSSION OF SELECTED PREDICTORS 

No attempt was made to justify rigorously, using either 
Meceorological principles or statistical testing, the analog 
elements that exhibited high indicated forecast potentials 
and thus were selected as predictors. There are, however, 
some obvious results and trends worth noting. See Appendix 
G for the selected predictors for each forecast period. 

Extrapolation of past intensity changes in operational 
forecasts has been used extensively and has been shown to 
produce reliable results for short-term forecasts. Of the 
derived elements, the one exhibiting the highest indicated 
forecast potential for the 24-hour forecast was the past 24- 
hour change in intensity. This element was also selected 
as a predictor for the 48-hour forecast period, but it ex- 
hibited a lower indicated forecast potential than it did 
for the 24-hour forecast period. For the 72-hour forecast 
period, no past intensity change elements exhibited large 
enough forecast potentials to warrant selection as predic- 
tors. This trend agrees with the fact that extrapolation of 
past intensity changes in operational forecasting has given 
more accurate results for the short-term forecasts. 

The results indicated that the height and latitude of 
the 700-mb subtropical ridge at the longitude of the cyclone 


Meee related to intensification during all forecast periods. 
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Of all the original and derived elements, these quantities 
and derivatives involving them appeared as predictors the 
most persistently. The resulting weight factors of predic- 
tors involving these quantities generally increased in the 
most extensive forecast period. For example, the weight 
factors derived for the 700-mb ridge latitude increased from 
1.12 and 1.13 for the 24- and 48-hour forecast periods to 
lee, for the 72-hour forecast period. The derived analog 
element involving the difference in 700-mb heights of the 
ridge and cyclone divided by the difference in latitude 
Showed a similar trend. It was not designated as a 24-hour 
predictor, but the weight factor associated with this ele- 
ment increased from 1.13 for the 48-hour period to 1.33 for 
the 72-hour period. The past 24-hour change in this pre- 
dictor also exhibited weight factors increasing with the 
forecast Dommeeece from tao EO 1,.37eto 1.81 in the 2?'r, 
48-, and 72-hour forecast periods, respectively. 

The past 12-hour speed of movement element displayed 
decreasing indicated forecast potentials from the 24- to 
ehour forecast periods and was defined as a predictor for 
only the 24- and 48-hour periods. The weight factors 
associated with this predictor decreased from 1.42 for the 
24-hour forecast to 1.07 for the 48-hour forecast. The 
past 24-hour speed of movement element was selected as a 
predictor for the 72-hour period, but possessed the lowest 


Beaght factor for that period. 
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mivemw,00-mb tmouch height sand the past 24-hour change in 
this quantity Beiaieaia the 48-hour forecast predictors, in 
that they possessed the two highest indicated forecast po- 
tentials and associated weight factors. The 700-mb trough 
height also possessed the highest indicated forecast poten- 
en of any of the original analog elements for the 72-hour 
forecast period. 

The estimated 1000/700-mb thickness element was not 
itself defined as a predictor for any of the forecast per- 
iods. However the past 24- and 48-hour changes of this 
ent did exhibit sufficient indicated forecast potentials 
to be designated as predictors. There was little observable 
trend in the indicated forecast potential of these predic- 
tors from one forecast period to the next. The past 48-hour 
change in thickness was defined as a predictor for the 24- 
hour period; no thickness elements were defined as predic- 
tors for the 48-hour period; the past 24-hour and again the 
past 48-hour changes in thickness were designated predictors 
for the 72-hour period. 

ie is bes interesting to note that the minimum /00-mb 
height above the cyclone was defined as a predictor for the 
24-hour forecast period only. However, derived elements 
involving this quantity were designated as predictors for 
all forecast periods. 

The elements designated as predictors from the results 


of this study are not necessarily considered to be the 
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mmmemal best predictors of tropical cyclone intensification 
among the possibilities derivable from the data. An almost 
fiini1te number of combinations of the original 18 analog 
elements could be evaluated as possible predictors. As will 
be discussed later, the results of just those elements con- 
sidered in this study could be altered by modified and im- 


proved procedures. 


Peer ORECAST VERIFICATION AND RESULTS 
1. Results of All Dependent Test Cases 

The analog forecasts derived for the 30 dependent 
and 61 independent test cyclones, as well as the available 
official JTWC forecasts supplied by FWC Guam were verified 
on the bastsmOtmmnvoumMain  eriteraas (1) absolute value of the 
Momeecast intensity errors; (2) the algebraic and absolute 
values of forecast intensity errors relative to the verify- 
ing intensities and expressed as percent errors. The first 
of these criteria is of little practical value in an opera- 
tional forecast procedure. It considers only the magnitude 
of the error without regard to the actual magnitude of the 
cyclone intensity. For the purposes of operational forecast- 
ing and contingency procedures of all concerned activities, 
a 20-knot error in a 24-hour forecast, for example, would 
be considered misleading and inadequate for a tropical 
storm with a verifying intensity of 40 or 50 knots. The 
Same error for the same forecast period, however, would be 


considered adequate for a super typhoon with a verifying 
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[imfeensity of 130 knots [3]. Therefore it is the relative 
merecast error that is of primary concern. The Joint Typhoon 
Warning Center has established relative error acceptability 
emerceria for different forecast time periods. The criteria 
for the 24-, 48- and 72-hour forecast periods are defined 
in [3] and are presented in Table I. The criteria are re- 
laxed for more extended forecast periods as contingency re- 
quirements are more flexible over ionger time periods [3]. 
Average values of the absolute intensity errors are provided 
only to indicate to what extent the error magnitudes vary 
from one forecast period to the next and to provide an addi- 
tional comparison oe analog forecasts to those of JTWC. 
inftemanalog forecast results are presented for all three 
forecast methods:modal intensification, intensification based 
on weighted predictors, and intensification based on weighted 
climatology (weighted combination of first two methods). 
Table II and Table III present the average absol te 
Forecast errors for all of the dependent and independent test 
cases, respectively. Generally for both set of cases, tie 
largest increase Ses from the 24- to 72-hour forecast 
periods occurred between the 24- and 48-hour forecast per- 
jods. Relatively little forecast deterioration was observed 
between the 48- and 72-hour periods. This trend was exhibi- 
ted by all forecast methods with two exceptions: the modal 
intensification for the dependent cases and the weighted 


climatology intensification for the independent cases. 


48 














Table IV and Table V provide the 24-hour relative 
forecast errors for the dependent and independent test 
cases. For the dependent cases the weighted climatology 
produced the best results, but for the independent cases, 
modal intensification was generally the most accurate. For 
the weighted predictors and weighted climatology intensi- 
aeons, the results deteriorated to varying degrees in 
gOing from the dependent to the independent test cases. 

But for the modal intensification, the results of the jnde- 
pendent cases showed improvement over those of the dependent 
cases, with the exception of the 10% error category. 

For both the dependent and independent test cases, 
the weighted predictors intensification showed the greatest 
Fraction of useful forecasts (< 30% error) that were accu- 
rate to the meaSurement of the Pere (< TOC error)’. For 
example, considering the dependent cases, (Table V), this 
Memecast method resulted in /54(50/67=.75) of one useful 
forecasts being accurate to the measurement of the data. 
These percentages for the modal and weighted-climatology 
Magensifications for the dependent cases were 55% and 622, 
respectively. This indicates that although the weighted- 
predictors intensification did not produce as many useful 
forecasts as the other methods, it became relatively more 
skillful than the other methods as the forecast error 
decreased. 

The 48-hour relative forecast errors are presented 


Mmretable VI and Table VIL. The modal intensification and 
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that of weighted climatology exhibited nearly the same re- 
eiets, especially for the useful (< 40% error) and adequate 
(< 30% error) forecast categories. All forecast methods 
suffered deterioration in going from the dependent to inde- 
pendent test cases. The modal intensification exhibited 
the largest fraction of useful forecasts that were accurate 
to the measurement of the data for both sets of test cases. 
The weighted predictors showed the smallest fraction for 
the dependent cases while the weighted climatology showed 
the smallest fraction for the independent cases. 

The relative forecast errors for the 72-hour fore- 
east period are presented in Table VIII and Table IX. For 
the dependent cases, the best results were observed for 
the modal and weighted climatology intensifications, while 
modal intensification dominated with the independent tests. 
Again the results generally deteriorated from the dependent 
to independent cases for all forecast methods, but to a 
dJesser extent for the modal intensification. The modal 
intensification again exhibited the largest fraction of 
useful forecasts hee were accurate to the measurement of 
the data for both sets of test cases. As with the 48-hour 
forecasts, the weighted predictors showed the smallest 
fraction for the dependent cases, and the weighted clima- 
tology showed the smallest fraction for the independent 
cases. The fraction associated with the weighted predictors 
showdd an increase from the dependent to independent test 


Mases (16/64 = .25 to 16/55 = .29). 
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Table X presents the relative forecast errors of the 
Ginmaependent test cases stratified by initial intensity. Two 
error categories are presented -- those that resulted in in- 
adequate forecasts, and those that resulted in forecasts 
accurate to the measurement of the data. The test cases 
were stratified on the basis of whether or not their initial 
intensities were of typhoon strength. For the 24- and 48- 
hour forecast periods the number of occurrences of inadequate 
forecasts was substantially less for the typhoon cases than 
it was for the storm or depression cases for each of the 
Maree FOrecast methods. But for the 72-hour forecast period 
the results were much more similar for the modal intensifi- 
cation, and essentially idanemee for the remaining two 
methods. The number of forecasts accurate to within data 
measurements for each forecast period was generally greater 
momeche typhoon casés than for the storm or depression 
Cases. For the 24-hour forecast period the number of for: - 
casts accurate to the measurement of the data was signifi 
cantly greater for the typhoon cases than for the depression 
Or eeorn fasesmror Each Lobrecast method. Except for the 
modal intensification, this stage-related trend decreased 
Substantially in the more extended forecast periods. These 
results imply that the intensifications of the more intense 
tropical cyclones are generally more predictable than those 
of the weaker, developing or dying, depressions and storms. 
Mr results in Table X are based on relative errors, and 


therefore do not imply that the absolute forecast errors 
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decrease with increasing initial intensity. But as stated 
previously, an absolute error of 20 knots is not nearly as 
significant for a 130-knot typhoon as it is for a 40- to 
50-knot storm. Therefore, if it is accepted that only rela- 
tive errors are meaningful, the implications of the results 
in Table X are valid. 

Frequency histograms of relative errors for the 
three forecast methods were constructed for both dependent 
and independent test cases. The histograms for the weighted- 
climatology forecast method are representative of all those 
constructed and are presented in Figures 1 through 6. In 
mmeecases, relative errors of + 54 or less were classified 
as zero forecast error. For positive errors, NZ error in- 
feeaedad all errors within (N-4)% and (N+5)% inclusive, and 
Morenegative errors, -N%Z error included all errors within 
Matt4)%, and (-N-5)Z, inclusive. Therefore the results 
presented in the eee do not exactly agree with those 
presented in the tables. 

The histograms for the dependent cases reveal that 
the Beecast ability of the weighted-climatology method 
decreased from the 24- to 48- to 72-hour forecast periods. 
That is, the numbers of occurrences of small relative fore- 
cast errors (< 20%) decreased and the distribution became 
more variable for the more extended forecast periods. No 
persistent skewness toward positive or negative errors 
resulted. However it was observed that the largest forecast 


errors were consistently positive errors. 
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The results for the independent test cases differed 
somewhat from the dependent cases. The most obvious differ- 
ence was the persistent skewness toward negative forecast 
Orrors, especially in the 24- and 48-hour forecast periods. 
The fact that this result did occur might have made possible 
the use of an adjustment factor or function to modify and 
improve the results. This was not attempted, however, in 
this study. The histograms reveal a decrease in the number 
of occurrences of small forecast error with increasing fore- 
cast periods. One similarity of these results to those of 
the dependent cases was the fact that the largest errors 


were consistently positive errors. 


Dis Comparison of Analog Forecasts to Official 
JTWC Forecasts 


As a further measure of the significance of the 
forecast results, comparisons of the analog forecasts to 
those of the official JTWC forecasts were made for homo- 
geneous sets of dependent and independent test cases. Com- 
parisons were made for the 24- and 48-hour forecast periods 
only due to the very limited number of applicable 72-hour 
forecasts available from FWC Guam. The analog results were 
recomputed based on only those test cases for which the 
official forecasts were available. 

The average absolute forecast errors are presented 
in Table XI and Table XII for the dependent and independent 
tests,respectively. As the figures indicate, the average 


errors resulting from all three analog methods were from 35% 
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Mmeeermthan, to at worst equal to those of the official fore- 
M@mets £oOr the dependent test cases. For the independent 
Mees, the errors for the analog forecasts were from 17% to 
mmost 60% greater than official forecast errors. The 
memented predictors method gave the worst results for the 
Imdependent cases. Neglecting this method, the largest per- 
cent increase in error was 33% (for the weighted climato- 
Meey). tthe overall best results, considering both forecast 
periods and both sets of test cases, were produced by the 
Modal intensification. However, the weighted predictors 
and weighted climatology did result in smaller average fore- 
cast errors for the 24-hour dependent forecasts than did 
the modal ap eERe eee 

Table XIII shows the comparisons of relative fore- 
cast errors among the analog methods and the JTWC forecasts 
for the dependent tests. These results reveal that all of 
the analog methods exhibited some equal or Superior foreca t 
mummetrcy Over the official forecasts. Exceptions were the 
modal and Pe. woe eer 24-hour forecasts for the 
< 30% error category, and the weighted-predictors 48-hour 
forecasts for the < 20% and < 10% error categories. All but 
O@imenos the analog methods exhibited a greatest fraction of 
useful forecasts that were accurate to the measurement of 
Mmemdata for both periods than did the official forecast. 
Mememindicates that aes the forecast error decreased, the 


analog methods (with one exception) became relatively more 
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feeebrul than did the official forecasts. The one exception 
to this result was the weighted predictors intensification 
for the 48-hour period only. However, the fraction exhibited 
by the weighted predictors method for the 24-hour forecast 
period was substantially greater than any of the others. 

The relative error comparisons for the independent 
Meee are shown in Table XIV. For the most part, the offi- 
cial forecasts proved to be superior to the analog fore- 
casts. Exceptions were the 24-hour weighted climatology 
Memeereasts for the < 10% category error category and the 
48-hour modal and weighted climatology forecasts for the 
Seo, e€rrer category. «For all analog methods for both fore- 
tee periods, the superiority of the official forecasts 
Gemereased with decreasing error categories from < 30% to 
Eee, errors. However, the trend is just the opposite in 
Pens from the < 40% to < 304 error categories. With one 
Pyeeeption, the fractions of - useful forecasts that were 
accurate to the measurement of the data were greater for 
the analog forecasts than for the official forecasts. 
Again, the Berea — the weighted predictor method for 


the 48-hour forecast ok-b ap Moke an 


C. CONCLUSIONS 
1. Implications of Forecast Results 
mie overall results of the forecast verification 
@re considered competitive with official forecasts and thus 


encouraging. The general deterioration of the results from 
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the dependent to independent test cases was expected; none- 
theless substantial numbers of useful and adequate analog 
forecasts were generated for these independent cases. The 
results were especially encouraging for the 48- and 72-hour 
forecasts, as these extended periods offer the forecaster 
the greatest challenge. If these results were maintained 
under conditions of operational forecasting, they would at 
least provide the forecaster with a starting point, other 
man linear extrapolation and climatology,for deriving his 
prediction. With a somewhat reliable starting point, the 
Panished forecast might well be more accurate. 

Of the three analog forecast methods applied, those 
involving the use of various analog elements possess the 
greatest potentials for improvement. Modal intensification 
might be improved slightly by experimenting with the screen- 
ing meieeria, but little else could be done to reduce the 
forecast errors. However, an almost unlimited amount of 
experimentation with the analog elements and their forecast 
potentials is possible. 

It has been established that the analog forecast 
methods performed more accurately on test cases used to 
derive the predictors and their associated weight factors. 
The analog elements associated with the small sample of 
test cases utilized for this purpose eee eedaey are not 
representative of all the element distributions. However, 


if all available past data were utilized in deriving the 








mmedtctors and their weight factors, the resulting distri- 
mmemens of the analogyeelements would bewvery representative 
of the possible values that each element could assume. 
Therefore, the likelihood would be high that all but extreme- 
ly uncommon cyclones (for which forecasts must be issued), 
would resemble at least some cyclones used to derive the 
predictors and weight factors. Consequently, the relation- 
ship of any given current cyclone to the information used 
in deriving the forecast scheme would approach the relation- 
ship between the dependent test cases and the information 
fereaein deriving their forecasts. Thus it is strongly 
believed that more accurate results could be obtained by 
increasing to the Serene extent the number of test cases 
merlized in deriving the predictors and their weight factors. 

As stated earlier, the number of elements derivable 
from the original 18 analog elements is enormous. Evalue- 
tion of elements other than those considered in this stuc 
might well lead to potentially better predictors and for: - 
east results. Also the technical procedures of element 
Me luation and forecasting could be modified, or completely 
different avenues of approaches followed. 

2. Evaluation of Investigative Procedures 

The encouraging results of the analog forecasts 
derived in this study strongly suggest that the methods 
Miamized in defining the predictors and computing associated 
weight factors possess at least some degree of skill. In 


Spite of some crude and intuitively derived procedures 
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involved in the method, the resulting predictors apparently 
do possess at least some degree of forecast potential. 
Examples of these crude procedures are the definitions of 
the acceptance intervals for original and derived element 
matching. A "sledge hammer” effort was also applied in the 
forecast method to alleviate the problems of variable 
chance element matching. 

The most severe shortcoming of the methods utilized 
was the failure to properly account for the variable proba- 
bility of chance element matching. The effects of this 
factor were not realized until forecast outputs were 
analyzed. However, it is likely that the same distorting 
effects were present in the element investigation procedures. 

In the second and third approaches to the main ob- 
jective, the selected analogs ee jase TH on the basis 
of a maximum allowable deviation in intensification from 
that of the test case. First, all analogs that intensified 
Meeenin +t 10 knots per 24-hour period were classified and 
Pie degree of Mant matching determined. Then all analogs 
Mee inderwent DeCHanze=riermecnsiiy of + 20 Knots, + 30 
Knots and finally + 40 knots per 24-hour period were classi- 
fied and for each maximum allowable intensification devia- 
tion, the degrees of element matching were computed. It is 
obvious that as these deviations increased, the number of 
acceptable analogs increased, which in turn, generally in- 


eCreased the number of possible analog element matches. 
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Ferma patrticular anmlog element, ae populasdoen could 
be defined as the entire set of values for that element con- 
Mained in the complete data file. The different classes of 
analogs represent samples of various sizes of this element 
population. Generally as the number of analogs increased, 
ieee increasing intensification range, the size of the ele- 
ment samples also increased. The variable sample sizes is 
Ome fLactor of the chance element matching. 

Tire other staceor 1s “tie location; within an element 
population, of the value of the test case analog element. 
Mmmact 1S, in general, the farther (closer) the test case 
memeue is from the expected (i.e. mean) value, the smaller 
Gereater) the probability of getting a match by pure chance. 

ine wre Pimoeceduresson this sttdy, the degree of 
Matching was defined as the percent of matches to the total 
number of possible matches. However, the ratio involved in 
ees computation can also be considered an average of ranks. 
memes ratio is mothing More than’ the sum of ranks (0 or 1) 
divided by the number of ranks. 

Tae Cemtral Limit PMhe@rem can be applived to this 
peoblem. The implications are that as the sample size in- 
creases, and as the test cyclone element value departs the 
expected value in the element population, the smaller the 
mecbability® that a specified average rank (e.g. .50) will 
Se@cur due to chance) @lone. This implies a smaller proba- 


Peerity ofWawspecificed degree of maitching (e.g. 50%) due to 
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chance as the sample size increases. But as the sample size 
increases and as the test case element value approaches the 
expected value in its distribution, the greater the proba- 
bility of a specified average rank due to chance. This 
implies a greater probability of a specified degree of match- 
ing,due to chance alone, as the sample size increases. The 
Significance of a match varies with the location of the test 
case element value in its distribution. Lo saeeount for this, 
the size of the acceptance interval should vary from smaller 
Magnitudes in high density regions of the population to 
larger magnitudes in low density regions. 

To properly pomenietatc for chance element matching, 
both the analog class sizes and the location of the parti- 
cular test case element value within its population must be 
accounted for. This requires a thorough study of the dis- 
mepoutions of all predictors to be utilized. 

Since chance matching could tend to increase the 
degree of matching for both large and small analog classes 
(depending on the location of the test case element value 
Mitts population), athe procedures taken Imes “Skidy £6 
@empensate for chance matching were only partially valid. 
The procedures were proper only when the test case element 
values were in low density regions of their distribution. 

Regardless of the shortcomings of the procedures in 
this study, the results warrant further study into the 
potential of analog schemes in forecasting tropical cyclone 


memeensification. Pies Lewewihat concentrated efforts in 








mies area could lead to improved and more significant re- 
sults, and the evolution of a reliable and efficient in- 
tensification forecast aid. 
Dee ihowghts on Subsequent Studies 

In all the investigative and forecast procedures 
ieea tized meni s study, the first operation was to screen 
out those six-hourly reports of past tropical cyclones that 
exhibited similarity to the test case with respect to loca- 
@men and time of occurrence. Further screening criteria 
then determined those analog candidates that resembled the 
test case with respect to the analog elements. The ideal 
analog procedure would be to continue the screening criteria 
until one analog candidate was defined to be the most simi- 
lar to the test case. This ideal procedure would require 
that the data observations be exact. It would also be 
necessary to know exactly what parameters are the ideally 
best predictors. These ideally best predictors would have 
to possess precise relationships to intensification, and 
the relationships would have to hold true for all cyclone 
irrences. In reality the data are not exact, and reliable 
predictors with precise relationships to intensification 
have not yet been defined. This forces the analog technique 
to rely on average behavior in order to generate a quanti- 
tative forecast. With the use of averages, it is hoped 
that random observational errors and non-persistent rela- 


tionships will be smoothed to a degree that enables the 
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utilization of the data in an adequate forecast scheme, 
Indeed, if the data were exact and precise relationships 
between predictors and intensification were known, deter- 
mMinistic techniques would replace statistical schemes. 

To generate a forecast, it would be necessary to simply 
substitute predictor values in an equation and solve for 
the intensification. 

It is felt, however, that some efforts should be 
made to strive for a more ideal analog forecast technique, 
involving more sophisticated screening criteria aimed at 
reducing the number of analog cyclones used in generating 
the forecast. A continuous~-type degree of matching 
criterion was attempted in the first approach of this study. 
Due to the inconclusive results of this procedure, a match, 
feematch criterion was investigated. A combination of 
these two matching criteria might lead to better results. 
That is for each analog element apply a match, no-match 
screening procedure to eliminate grossly unmatched values, 
and to identify those values within some prescribed range 
from the test case ee. With a large degree of noise 
eliminated, then proceed with a continuous-type degree of 
matching operation. For example, the difference in element 
values relative to the magnitude of the acceptance interval 
(used in the first step) could be examined. This type of 
procedure could possibly define more precise and persistent 
relationships between analog elements and intensification 


than those which were defined in this study. With the use 
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of more precise relationships, the screening procedure in 
the forecast scheme could be more restrictive and thereby 
select a fewer number of analogs that are more similar to 
the test case with respect to predictor vajues. In any 
event, further studies might reveal to what extent average 
behavior must be relied upon to smooth out the data errors 


and non-precise and quasi-persistent relationships. 








TABLE. I 


Forecast-intensity relative error acceptability criteria, 
as defined by JTWC. 


E = forecast-intensity error relative 
EG) Vermirteqtion intensity 





Forecast Value Forecast Period 


24-hour 48- and 72-hour 








Accurate to within 
intensity measure- 


Ment error Bos 02 ES E04 
Adequate Ee 07 : oe Ser 
Useful E < 30% E < 402 
Inadequate E307 Ee 407 














TABLE Il 


Pwemape absolute forecast-intensity errors (kt) for all de- 
pendent test cyclones. Numbers in parentheses indicate the 
number of applicable forecasts for each forecast period. 





Forecast Forecast Method 
Period . 
(Hr) Modal Weighted Weighted 
Predictors Climatology 

24 (30) 16 de5 13 

48 (29) 16 23 20 

aZ (25) 19 23 Dib 

TABLE III 


Average absolute forecast-intensity errors (kt) for all in- 
dependent test cyclones. Numbers in parentheses indicate 
the number of applicable forecasts for each forecast period. 








Forecast Forecast Method 
Period . 
Cia) Modal Weighted Weighted 
Predictors Climatology 
24 (61) 16 nS 26 
48 (59) 26 28 28 
woe C5 1) 28 5 32 














TABLE IV 


Percent of all dependent test cyclones exhibiting the indi- 
mated relative error (E) in the 24-hour forecast intensity. 
Results based on 30 forecasts. 


E = forecast-intensity error relative to verification 


hie ens ity. 


Forecast 
Method - E(Z) 


Modal | 87 73 57 40 

Weighted 

Predictors 83 67 62 50 

Weighted ; 

Climatology 90 80 70 50 
TABLE V 


Percent of all independent test cyclones exhibiting the 
[Mmaicated relative errors (E) in the 24-hour forecast 
intensity. Results based on 61 forecasts. 


EB = forecast-intensity error relative to verif-cation 
Mircensity. 


Forecast 
Method E(Z) 








Modal 90 82 64 sul 
Weighted 
Predictors 79 66 54 38 
Weighted 
Climatology 84 74 54 36 











ib LE ev TL 


Pemecnt of all dependent tests eyclones exhibiting the indi- 

mmmed relative errors (E) in the 48-hour forecast intensity. 

Reisults based on 29 forecasts. 

Mmeeeorecast-intensity error relative to verification 
Mmeensity. 


herceast 
Method E (2) 


Modal 79 70 59 41 
Weighted 
Predictors 70 59 43 zal 
Weighted 
Climatology 76 70 59 34 


TABLE VIL 


Percent of all independent test cyclones exhibiting the 
[mmemeated relative errors (E) in the 48-hour forecast 


intensity. Results based on 59 forecasts. 
B= forecast-intensity error relative to verification 
matensity. 





Forecast 
Method E(Z) 








Modal TAL Sa 44 23 
Weighted 
Predictors 64 49 30 20 
Weighted 
Climatology 68 54 34 15 











TABLE VIII 


Percent of dependent test cyclones exhibiting the indicated 
relative errors (E) in the 72-hour forecast intensity. 
Results based on 25 forecasts. 


P= tarecast-intensity error relative to verification 
imeensity . 





Forecast 
Method ECA) 


Modal 68 56 5 40 

Weighted 

Predictors 64 56 39 16 

Weighted 

Climatology 80 56 36 24 
TABLE Lx 


Percent of all independent test cyclones exhibiting the 
mearcated relative errors (E) in the 72-hour forecast 
intensity. Results based on 51 forecasts. 


E = forecast-intensity error relative to verification 
fmtensity. 


Forecast 
Method E(Z) 





Modal 67 oS) 43 Zee 
Weighted 
Predictors 315: 41 a2 16 
Weighted 
Clanat ology 59 a5. eal 16 
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Figure 1. Frequency historgram of 24-hour relative fore- 
cast errors for weighted climatology forecast method, 
dependent test cyclones. 
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Figure 2. Frequency histogram of 48-hour relative forecast 
errors for weighted climatology forecast method, dependent 
test cyclones. 
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Figure 3. Frequency histogram of 72-hour relative forecast 


errors for weighted climatology forecast method, 


test cyclones. 
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Figure 4. Frequency histogram of 24-hour relative forecast 
errors for weighted climatology forecast method, independent 
test cyclones. 
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Figure 5. Frequency histogram of 48-hour relative forecast 


errors for weighted climatology forecast method, independent 
test cyclones. 
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Figure 6. Frequency histogram 
errors for weighted climatology forecast method, independent 


test cyclones. 
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APPENDIX A 
Original Analog Parameters and Error 
Estimates of the Analog Elements 

The following Pe omere ys were compiled at six-hourly 
intervals for use with the TYFOON analog forecast scheme. 
Estimated maximum observational errors in the analog 
elements were provided through informal discussion with 
Commander J. D. Jarrell, USN. Commander Jarrell is a co- 
orjginator of the TYFOON analog forecast scheme and has 
served as Operations Officer for the Joint Typhoon Warning 
Center, Guam. Error estimates for the past movement ele- 
ments are based on estimated errors of six nautical miles 
in the latitude and longitude of a cyclone position. 
Identification Elements 

ili. Storm name (a number is given unnamed storms) 


Zoe stOrmerdentity number, Month, year and cyclone 
sequence number 


3. Observation number 
oo) Lear 

Bie Month/Day 

6. Hour (GMT) 

7. \Latitude 


8. Longitude 
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Analog Elements 


1. 


10. 


i. 


2 . 


14. 


15 


16. 


Estimated Maximum 


Post, L2Z=“hr direct1 on 
of movement 6 


Past 12-hr speed 
of movement 1 


Past lG=nr Gireceion 
of movement < 6 


Past 18-hr speed of 
movement <1 


Pase —24—-hr direction 
of movement <0 


Past 24-hr speed of 
movement << l 


Past 48-hr direction 
of movement < 6 


Past 48-hr direction of 
movement , < 1] 


Size (average radius of 
outer closed isobar in 


whole degrees of latitude) 2 


Past 12-hr change of size 
in degrees latitude 4 


Minimum observed sea-level 
pressure in whole millibars 5 


Past 12-hr change in minimum 


Observatlonal Brrors 


degrees 


Kiet 


degrees 


knot 


degrees 


Knot 


degrees 


knot 


degrees latitude 


degrees latitud 


Magill pay s 


sea-level pressure 10 millibars 
Maximum wind 10%Z of the observed 
value 


Minimum 700-mb height 


(above the cyclone) 10 meters 


Latitude of the 700-mb 
ridge north of the cyclone 2 


700-mb height at the ridge 
line north of the cyclone Z0 
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degrees latitude 


meters 





7. 


18. 


Pongutude at SON Jeo beeine 
nearest /00-mb trough 
west of the cyclone 


700-mb height at the 


imtersection of the trough 
line and 35N 
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2 degrees longitude 


30 meters 





AMeEND IX 2 


Dependent Test Cyclones 


Dependent test cyclones. The following tropical cyclone 
reports were selected to examine the forecast potential of 


the analog elements. 


Name MO-DA-YR Time initial 
(GMT) Intensity 
(kt) 
Betty 8-06-69 1200 45 
Gora 8-17-69 0000 50 
Doris 8-02-69 0600 30 
Wendy 8-29-68 1800 100 
Wendy 8-31-68 0000 140 
Shirley 8-17-68 1800 30 
Shirley 8~18-68 1200 50 
Agnes 8-30-68 0600 45 
Agnes 9-01-68 0000 85 
Agnes 9-03-68 0600 130 
Agnes 9-04-68 0600 145 
Agnes 9-05-68 0600 aes: 
Agnes 9-07-68 0600 70 
Agnes 9-08-68 0600 60 
Polly 8-07-68 1200 40 
Polly 8-09-68 0000 55 
Polly 8-11-68 0000 50 
mod ly 8-13-68 - 0000 Sie 
Louise 8-20-67 0000 30 
Louise 8-21-67 0000 55 
Marge 8-26-67 1200 100 
Opal © 9-01-67 0600 we) 
Opal 9-03-67 1200 155 
Opal 9-11-67 1200 90 
Rita 8-04-66 1200 60 
Rita 8-05-66 1200 15 
Carmen 8-10-63 0000 40 
Carmen 8-12-63 1800 265) 
Carmen 8-13-63 0600 110 
Carmen 8-14-63 0000 85 
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APPENDIX C 
Acceptance Intervals Defined for the 
Original Analog Elements 
The acceptance intervals defined for the 18 original 
analog elements are listed below. The range of each in- 
terval is defined as plus or minus one-fourth of one 
standard deviation of the respective element's distribu- 
tion over the entire data file obtained for the study. 
This range is centered on the test cyclone element value. 
If this fraction of the standard deviation resulted in a 
number less than one, the acceptance interval was defined 
wet. 0. Intervals vice discrete values were preferred 
in order to prevent the necessity of identical element 


Mewes in order for a match to occur. Refer to Appendix 


A for analog element definition. 


Analog Element Number Acceptance Interval 
1 i? seedies recs 
2 + 1.4 knots 
3 te) Omdiee Tees 
4 Jo oO. S 
5 + 7.6 degrees 
6 tee ono CS 
7 + 7.3 degrees 
8 mp ho 8 henleie cs 
9 ft te@edesmees, latitude 
10 fo mdernrces latitude 
iLek a2. / Molde bars 
eZ to) 2 eee Ee bait S 
iS + 7.9 knots 
14 + 4.7 dekameters 
15 + 1.4 degrees latitude 
16 + 1.0 dekameters 
IY Peoemeaezrees longitude 
18 + 1.0 dekameters 
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APPENDIX D 
Indicated Forecast Potentials of the 
Original Analog Elements 

Indicated forecast potentials of the original analog 
elements for each forecast period are listed below. The 
indicated forecast potentials are expressed as a percent 
Seeeene total number of possible occurrences in which the 
trend of decreasing degree of matching with increasing 
forecast error occurred for each element. The results are 
based on thirty 24-hour forecasts, twenty-nine 48-hour fore- 
casts and twenty-five /2-hour forecasts. Refer to Appendix 


meror element definitions. 


Analog Element Forecast yheotential for 
Number Given ‘Forecast Periods (%) 
24-hr 48-hr 72-hr 
1 nny Zo 8 
Z 41 a2 12 
3 ily 29 eZ 
4 MeN 18 8 
5 25 30 “al 
6 2a 19 ig 
i 24 29 14 
8 Ly iy 2 
9 She) 28 BZ 
10 3 Mais 28 
le 25 31 16 
eZ 47 Ze 24 
: 3 ay te 20 
14 53 al 20 
15 37 34 ie 
16 40 38 S\H4 
17 \S) 28 20 
18 20 41 36 
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APPENDIX E 
Derived Analog Elements and their 
Acceptance Intervals 
The derived analog elements, formed by various combi- 
nations of some of the original analog elements, are listed 
below together with the acceptance interval that was 
selected for each. The Pre: interval fer each ele- 


ment is centered on the test cyclone value of that element. 


Derived Analog Acceptance 
Elements Intervals 
1. (700-mb ridge height-mininum ieee o To 
700-mb height) 
2. (700-mb ridge lat-cyclone lat) + 2 deg lat 


3 mOO-mb ridge height-minimum 700-mb heteise 
700-mb ridge lat-cyclone lat 


+ 20 m/deg lat 


4. (700-mb trough height-minimum 


700-mb height) 20 m 


fete 


5. (700-mb trough long-cyclone long) + 2 deg long 


6 700-mb trough height-minimum 700-mb height 
° ~ 700-mb trough long-cyclone long 


+ 10 m/deg long 
7. 1000/700-mb thickness + 2 deg lat 


S$. Past 24-hr change in derived 
element #2 toe2gedeg lat 


9. Past 24-hr change in derived 


éhement #3 10 m/deg lat 


Near 
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Oe 


re 


2. 


13. 


14. 


ieee 


16. 


17. 


iro 


ieee 


20. 


2 ANE 


15 5 


ZS) 


24. 


25. 


26. 


21. 


Past 24-hr change 
element #5 


Past 24-hr change 
element #6 


Past 24-hr change 
element #7 


Past 24-hr change 
(maximum wind) 


Past 24-hr change 
pressure 


Past 24-hr change 
700-mb height 


Past 24-hr change 
Brough height 


Past 24-hr change 
Yidge height 


Past 48-hr change 
element #2 


Past 48-hr 
element #3 


change 


Past 48-hr 
element #5 


change 


Past 48-hr change 
element #6 

Past 48-hr 
element #7 


change 


Past 48-hr change 
(maximum wind) 


Past 48-hr change 
pressure 


Past 48-hr change 
700-mb height 


Past 48-hr change 
Breoue hn height 


Past 48-hr change 
ridge height 


in 


ante 


a0 


in 


in 


aL ea 


aL 36) 


in 


in 


Ty 


a1 Fi 


on 


in 


aL) 


iLie" 


in 


in 


derived 


derived 


derived 


intensity 


sea-level 


minimum 


70O0—mb 


700-mb 


derived 


derived 


derived 


derived 


derived 


intensity 


sea-level 


minimum 


700-mb 


700-mb 
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10 


10 


10 


10 


10 


HO 


10 


10 


10 


10 


10 


10 


JELe 


deg long 


m/deg long 


kt 


mb 


deg lat 


m/deg lat 


deg long 


m/deg long 


kt 


mb 





APPENDIX F 
Indicated Forecast Potentials of the 
Derived Analog Elements 

Indicated forecast potentials of the derived analog 
feements for each forecast period are listed below. fMIhe 
indicated potentials are expressed as a percent of the 
total number of possible occurrences in which the trend 
of decreasing degree of matching with increasing forecast 
Seror occurred for each element. Refer to Appendix E 


mes element definitions. 


Analog Element Forecast Potential for 
Number Given Forecast Periods (%Z) 
24-hr 48-hr 72-hr 

1 29) : a1 16 

2 20 34 20 

3 23 34 28 

4 20 a4 16 

5 exe) 31 20 

6 40 aie 23 

7 27 24 Peer 

8 29 26 25 

9 ; 39 41 38 

10 eel. 33 21 

ll 18 26 17 

12 29 2 25 

13 46 37 ZA. 

14 32 Bee 2 1. 

15 18 ay, 13 

16 29 48 21 

tf 36 a3 2 Al 

18 a7 Wy 14 

19 36 25 18) 

20 28 29 18 

21 20 dey. 9 

22 36 2 Pad 

23 20 2 18 

24 28 33 18 

75 8 4 0 

26 24 29 23 

7 | 24 29 9 
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APPENDIX G 


Selected Predtetors=tor cach torecast Period 


The original and derived analog elements that were 
selected as predictors for each forecast period are listed 


below along with their associated weight factors. 


24-hr Forecast Period 


Preaictor Wetehteractor 
Past 12-hr change in sea-level pressure few 2 
Past 24-hr change in intensity (maximum wind) hee SY 
Past 12-hr speed of movement de 4 

700-mb trough height-minimum 700-mb height 1.21 

700-mb trough long-cyclone long ° 
700-mb ridge height diez 
Past 24-hr change in 

poO0-mb ridge height-minimum /00-mb height 1.18 

700-mb ridge lat-cyclone lat : 
700-mb ridge lat eels? 
Past 24-hr change in 700-mb ridge height 1.09 
Past 48-hr change in thickness 1.09 
Past 48-hr change in 

700-mb ridge height-minimum /00-mb height 1.09 

700-mb ridge lat-cyclone lat : 
700-mb trough long i..0'0 
minimum 700-mb height 1.00 


size i0'0 


48-hr Forecast Period 


Beedictor Wememt Factor 
Past 24-hr change in trough height 1.60 
Past 24-hr change in 

move-mb ridge height-minimum /00-mb height 1.37 

700-mb ridge lat~-cyclone lat ° 
700-mb trough height eo 
700-mb ridge height oe 
current intensity (maximum wind) le 27 
Past 24-hr change in intensity (maximum wind) ee 3 
Past 24-hr change in minimum 700-mb height eS 
700-mb ridge lat 5 ee bs 
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Go-hWebeLlecact Period (continued ) 


meedictor Weenie ae © Or 
700-mb ridge height-minimum 700-mb height 1.13 
700-mb ridge lat-cyclone lat 
mood-mb trough height-minimum /00-mb height 1.13 


700-mb trough long-cyclone long 
Past 24-hr change in 700-mb trough 


long-cyclone long 1.10 
Past 24-hr change in 700-mb ridge height 0 
Past 48-hr change in sea~level pressure 0 
Past 12-hr change in speed of movement 1,07 
Past 24-hr change in direction of movement ie te 


¥2-nr Forecast Period 


Peedictor We tees Pact or 
Past 24-hr change in 
700-mb ridge height-minimum 700-mb height 1.81 
700-mb ridge lat-cyclone lat ; 
700-mb trough height ; ay es 
700-mb ridge lat eZ 
Size ee) 2 
Past 12-hr change in size TRS es 
muo-mb ridge height-minimum 700-mb height 1.33 
700-mb ridge lat-cyclone lat ; 
Past 48-hr change in thickness 6 8) 
Past 24-hr change in thickness 19 
Past 12-hr change in sea-level pressure 1.14 
Past 24-hr speed of movement 1.00 
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APY ENDS i 


Independent Test Cyclones 


Independent test cyclones. The 


following tropical 


cyclone reports were selected as independent test cases 


for the analog forecast progran. 


Name 


Cora 
Cora 
Pouly 
DEL x 
Trix 
eae xX 
Bess 
Bess 
Wendy 
Wendy 
Wendy 
Hope 
Hope 
Marge 
Marge 
Opal 
Opal 
Opal 
Rita 
Rita 
Rita 
Alice 
Alice 
Alice 
Alice 
Alice 
Tess 
Tess 
Tess 
Cora 
Cora 
Cora 
Cora 
Cora 
Kim 
Kim 
Kim 


MO pA— yk 


8-18-69 
8-20-69 
8-04-68 
8-24-68 


8-25-68 - 


8-26-68 
9-01-68 
9-03-68 
9-02-68 
9-04-68 
9-06-68 
8-05-67 
8-06-67 
8-24-67 
8-25-67 
9-05-67 
9-06-67 
9-07-67 
8-02-66 
8-03-66 
8-06-66 
8-27-66 
8-29-66 
8-31-66 
9-01-66 
9-02-66 
8-12-66 
8-13-66 
8-14-66 
8-31-66 
9-01-66 
9-02-66 
9-04-66 
9-05-66 
8-04-65 
8-05-65 
8-06-65 


Gal 


Time ita tka 1. 
(GMT) Intensity 
(kt) 
1200 65 
0600 80 
0600 315) 
0600 35 
1200 40 
0600 50 
0000 40 
0600 50 
0000 OS 
0000 85 
0600 60 
0600 45 
0000 55 
0600 30 
0000 40 
0600 110 
0600 115 
0600 o5 
1200 25 
1200 35 
0600 50 
0000 5) 
1800 75 
1800 100 
1200 130 
1200 130 
UOMO. 25 
1200 35 
1800 25) 
0600 50 
1200 100 
1800 130 
0000 130 
0600 140 
0600 AS 
0600 50 
0600 60 





Name 


Lucy 

Lucy 

Lucy 

Lucy 

Lucy 

Olive 
Olive 
Olive 
Kathy 
Kathy 
Kathy 
Kathy 
Kathy 
Kathy 
Kathy 
Kathy 
Kathy 
Marie 


- Marie 


Marie 
Sally 
Sally 
Sally 
Sally 


MO-DA-YR 


8-16-65 
o=1 o 
S=i 7c 
o-20>o5 
Se 2 eo 
8-29-65 
8-30-65 
8-30-65 
8-12-64 
8-13-64 
8-14-64 
8-15-64 
Oyo 
8-18-64 
8-19-64 
8-20-64 
8-22-64 
O- 15-64 
8-15-64 
8-16-64 
9-04-64 
9-05-64 
9-06-64 
9-07 64 


She 


Time Meee te a 


(GMT) Intensity 
(kt) 
0000 85 
0000 2 O 
0000 120 
0600 EO 
0600 120 
0000 10 
0000 130 
1800 1G 
0600 35 
0600 80 
1200 80 
1800 60 
0000 7s. 
0600 110 
1200 138) 
1800 100 
0000 85 
0600 30 
1800 30 
0600 30 
0000 55 
0000 85 
0600 13:0 
* 1800 160 
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